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A purely ionic interpretation of the tilting mode insta- 
bilities in La2-2,Aa;Cu04 (A=Sr,Ba) is shown to be not self- 
consistent: the dominant factor influencing the doping de- 
pendence of the interlayer mismatch is the large change in 
the Cu-O bond length. But such a large dependence of bond 
length on valence leads to a strong electron-phonon coupling, 
contradicting the assumption that covalent effects were neg- 
ligible. This coupling is closely related to the vHs-JT effect. 
This new insight clarifles the role of the tilt-mode instabil- 
ities. The main JT coupling is not to these modes, but to 
the in-plane 0-0 bond stretching modes which split the vHs 
degeneracy. However, as these modes soften, they couple to 
the lower-lying tilt modes, so that the ultimate instability has 
a finite tilt component. 

This in turn simplifies the description of the JT instabili- 
ties. The tilt modes, while strongly coupled to the electrons, 
are quadratic in the tilt angle. On the other hand, the bond 
stretch modes have a large, linear coupling to electrons, with 
clear polaronic effects. This is connected with the fact that, 
near a vHs, valence fluctuations are slow and relatively long 
lived, coupling to local phonon modes. A striking result of 
this is that there will be a large polaronic band narrowing 
near the vHs, whether or not the vHs is near the Fermi level. 
This vHs-localized band narrowing provides a natural expla- 
nation for the common occurence of extended vHs's. 



I. INTRODUCTION 

Recent photoemission studies have demonstrated that 
a van Hove singularity (vHs) falls very close to the Fermi 
level in optimally doped BiaSraCaCuaOg (Bi-2212) |] 
and YBa2Cu307 (YBCO) [|. This lends strong support 
to the many proposals that the vHs plays an important 
role in high-Tc superconductivity |^ in these materials. 
The vHs's also couple to many of the phonon modes in- 
volved in the various structural phase transitions found 
in La2-a;Sra;Cu04 (LSCO) and related compounds 
which leads to a competition between superconductivity 
and structural transitions in the cuprates. 

However, these structural phase transitions are com- 
mon to most perovskites, many of which are insulating, so 
the notion that electron-phonon coupling plays an impor- 
tant role has been questionned. The conventional view 
of these transitions has been that they are purely ionic, 
driven by the Madelung energy or interlayer strain mis- 
match, and have little or nothing to do with the electronic 
system. 



This conventional picture has been modified by recent 
analyses. The perovskites have two classes of distortions, 
ferroelectric (FE), which generate a dipole moment in the 
unit cell, and antiferrodistortive (AFD), involving tilting 
of local octahedra. Cohen Q has shown that there can be 
large covalent contributions to the FE instabilities, even 
when the materials are insulating, if the Fermi level lies 
in a hybridization gap between two atomic species. (A 
similar result is discussed in Appendix A4 of Ref. 
This has stimulated a great deal of research showing the 
important role of correlation effects |^ in FE instabili- 
ties. Of particular interest is a recent two-band Hubbard 
model of the FE instability, introduced by Egami, et al. 
|lO| , which involves a one-dimensional (ID) version of 
the three-band Hubbard model regularly applied to the 
cuprates. 

The structural phase transitions in the cuprates are 
generally AFD, involving tilting of the CuOg octahe- 
dra. Such AFD instabilities are common in other per- 
ovskites, and indeed often compete with FE instabili- 
ties jl^. However, these AFD instabilities are generally 
found to be ionic rather than covalent - the tendency 
towards instability scales linearly with a tolerance fac- 
tor, t, which measures the degree of ionic size mismatch 
(]lT| (see, e.g., Eq. 1, below). Hence, the present paper 
provides a reassessment of the role of electron-phonon 
coupling in the structural instabilities of the cuprates. 

The paper is organized around three calculations. 
First, a mixed ionic-covalent model of the structural tran- 
sitions is recalled. It is shown that the same model 
which explained the T-dependence of the low temper- 
ature orthorhombic (LTO) phase in LSCO 1^ can also 
explain its doping dependence (Section II and Appendix 
A). [Strictly speaking, the dominant instability is to a 
phase of short-range order - identified with the pseu- 
dogap phase - but the transition temperatures for this 
phase and the LTO phase are proportional to each other.] 
While the model must contain a significant ionic contri- 
bution, there can also be a significant covalent contribu- 
tion. Secondly, the purely ionic model for the cuprates 
is analyzed in more detail, and a number of discrepan- 
cies are pointed out. Most importantly, it is found that 
the dominant factor comtroUing the doping dependence 
of the LTO transition is the large compressibility of the 
Cu02 planes: in ionic terms, there is a large decrease in 
the Cu-O distance when O^^ . This analysis is 

presented in Section HI and Appendix B. 

This large change of interionic spacing with valence is 
a signature of large electron-phonon coupling, as found in 
rare earth compounds p^ . Hence, a purely ionic model 
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with no electron-phonon coupling is not self-consistent. 
The contraction of the 0-radius couples to a number of 
phonon modes, but predominantly to in-plane modes, 
and not to tilting modes. This leads to a reanalysis 
of earlier results. It is now suggested that the major 
phonon coupling is indeed with planar modes, which cou- 
ple quadratically to the much lower frequency octahedral 
tilt modes. This weak coupling ensures that when the 
highly coupled bond stretching modes start to soften, the 
tilt modes go unstable first. 

This has important consequences for the vHs-JT ef- 
fect. The electron-tilt phonon coupling is quadratic in 
the tilt angle, and hence in the phonon creation operator 
- a strong violation of Migdal's theorem. In contrast, the 
the 0-0 bond stretching mode has a linear JT coupling 
to the vHs, via the length dependence of the hopping pa- 
rameter, tcuO- The soft mode behavior of, in particular, 
the 0-0 stretch modes is discussed in Section IV and 
Appendix C. Moreover, there is a clear polaronic contri- 
bution: as the hole hops around, it carries this lattice 
distortion with it. These polaronic effects are discussed 
in Section V and Appendix D. By estimating the mag- 
nitude of the polaronic corrections to the bandstructure, 
it will be possible to make connection with various pola- 
ronic or bipolaronic theories of superconductivity. 

One striking result is that polaronic effects can explain 
the existence of extended vHs's, even when the vHs is 
away from the Fermi level. Band narrowing effects are 
localized near a vHs due to the fact that, near a vHs, 
the material has an ioniccovalent crossover. Electronic 
motion is slow compared to phonon motion, so valence 
fluctuations are well defined (Section VI). Section VII 
summarizes some conclusions, while Appendix E applies 
some of these results to spin-orbit coupling. 

Once it is recognized that strong electron-phonon cou- 
pling effects can play an important role in the structural 
instabilities in the insulating perovskites, it then becomes 
plausible to assume that they can also drive a competi- 
tion between structural and superconducting instabili- 
ties, when the perovskites are doped. Indeed, undoped 
SrTiOs has an antiferrodistortive instability associated 
with the tilting of the TiOe octahedra, which is closely 
analogous to the LTO and low temperature tetragonal 
(LTT) phases of LSCO, whereas doped SrTiOs is super- 
conducting - one of the first known members of the fam- 
ily of high-Tc perovskites. A similar competition between 
tilt-mode instabilities and superconductivity is found in 
BaKi_3;Bi2,03 [|l3|. Interestingly, this is accompanied by 
softening of the ID 0-0 bond stretching modes, which 
also play an important role in the cuprates, as will be 
demonstrated below. 



II. STRUCTURAL PHASE TRANSITIONS 



A. Phase Transitions and Pseudogaps 

In LSCO/LBCO and related compounds, a number 
of structural phase transitions are observed. From high 
temperature to low, the principal phases observed are 
a high-temperature tetragonal (HTT) phase, where the 
octahedra appear to be untilted, a low-temperature or- 
thorhombic (LTO) phase, in which the average tilt is 
about an axis at 45° to the Cu-0 bond, and a low- 
temperature tetragonal (LTT) phase, where the tilt axis 
is along the Cu-0 bond. By substituting Nd for La, it 
is possible to generate an intermediate (Pccn) phase in 
which the tilt axis is at an angle between and 45° to 
the Cu-0 bond. The transition temperatures are a strong 
function of doping, with the LTO transition temperature 
monotonically decreasing with hole doping, while a full 
LTT transition is found only in LBCO near a doping 
x=0.125, or in Nd substituted LSCO. 

I have suggested that these transitions are represen- 
tative of a dynamic Jahn- Teller (JT) transition. The 
LTT phase would involve a predominantly static distor- 
tion, the LTO a dynamic tunneling between two adja- 
cent LTT tilts, and the HTT tunneling among all four 
tilts. The Pccn phase could be explained as an inter- 
mediate phase, where there is some LTO tunneling, but 
still a static LTT component. The HTT^LTO tran- 
sition is a transition between two dynamic JT phases, 
driven by intercell coupling associated with corner shar- 
ing of planar O's and long range strain forces. If the 
microscopic average local tilt is unchanged between the 
LTO and HTO phases, there may be little change in the 
electronic properties at the transition - this appears to 
be the case experimentally. However, it is possible that 
there is an additional transition within the HTT phase, 
associated with the onset of a tilting instability, which 
does couple more strongly to electronic degrees of free- 
dom. Such a crossover from the tunneling phase to an 
undistorted phase is a common feature of many dynamic 
JT transitions . 

This dynamic tilting onset can be identified with the 
'pseudogap' phase recently found by Hwang, et al. |l^] 
in LSCO, at a temperature T*, Fig. 1. Below this tem- 
perature, there is evidence for significant changes in the 
electronic properties of the material: the susceptibility 
starts to drop below T* , and the anomalous temperature- 
dependence of the Hall effect first appears. I have sug- 
gested that this anomalous temperature dependence fol- 
lows from the gradual opening of the short-range-ordered 
dynamic- JT phase structural gap, and signals the grad- 
ual crossover from the large Fermi surfaces to the small 
pockets seen in transport measurements. Hwang, et al., 
found a striking scaling of the pseudogap properties - 
the curves at different dopings all scaled to a universal 
behavior when plotted as functions of T/T*. The LTO 
phase transition has a similar scaling, with T^to — T* /2 
for all dopings. 

The dynamic JT model is consistent with frozen 
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phonon calculations by Cohen, et al., which find that in 
LBCO the HTT phase is unstable, and the LTO phase 
metastable with respect to the LTT phase A simi- 
lar pseudogap is found in YBCO, associated with peaks 
in the susceptibility and heat capacity jl^. This is 
associated with the onset of a pyramidal tilting instabil- 
ity seen in neutron diffraction [l9| ] and related phonon 
anomalie s [po| - |22| . This was originally interpreted as a 
spin gap ]23[|, possibly with spin-phonon coupling , 
but the recent evidence is suggestive of a gap in both 
charge and spin channels. 

Hence, I would identify this pseudogap phase with the 
onset of the JT-tilting instability. Thus, the calculations 
of this paper should refer to the T* transition. Since, 
however, most of the experiments are on the LTO phase, 
I will make use of the scaling found by Huang, et al. to 
compare the calculations for T* with experimental data 
for Tlto ^ T*/2. 

It should not be overlooked that such a dynamic JT 
effect is itself a potentially important (and highly anoma- 
lous) source of low energy excitations, which could act as 
additional bosons for causing electrons to pair up. As 
early as 1941, Teller was speculating that the dynamic 
JT effect could provide the driving mechanism for (con- 
ventional) superconductivity p5| . 

B. Doping Dependence of LTO Transition 

In the ionic picture for the cuprates, the structural in- 
stabilities are driven by interlayer bond-length mismatch 
p6| . The unit cell in the Cu02 layer is larger than that 
in the LaO layer, with the degree of mismatch measured 
by the tolerance factor 

t ^ dLao/V2dcu02: (1) 

where dLaO is the effective LaO distance in the LaO 
plane, with a similar definition for dcu02 ■ Commensura- 
bility strains will arise in LSCO whenever t is different 
from unity. Mismatch problems become more severe at 
lower temperatures, gradually causing the Cu02 layers to 
buckle - either by octahedral tilts in LSCO or by layer 
dimpling in YBCO (there may also be a pyramidal tilt in- 
stability in YBCO). Indeed, the structural instability can 
be interpreted in a purely ionic model, with the transition 
temperature proportional to the degree of mismatch. 

In undoped La2Cu04, the LTO transition can equally 
well be described by either a purely ionic model or a 
model with considerable electron-phonon coupling |^] 
(due to the strong correlation effects, this would be a 
spin-Peierls coupling, although the distinction is not ob- 
vious at the mean-field level). However, the transition 
could not be interpreted purely in terms of Fermi surface 
nesting; the 'bare' phonon frequency wq was always found 
to be imaginary, stabilized by anharmonic effects. This 
imaginary frequency could be caused by a combination 



of covalency effects and interlayer strains, 

Lo^„=ul-anh-T*\T+\ (2) 

(see Eq. 32c of Ref. Q ) . Here Uh = 1 + x is the number 
of holes, r+ is the interlayer strain, proportional to the 
lattice mismatch, is the bare frequency, and F* and a 
are appropriate coupling constants. The middle term of 
Eq. 2 is a covalency effect. Note that there is no electron- 
hole symmetry within the antibonding band: the maxi- 
mum frequency reduction corresponds to an empty band, 
Uh — 2. The electron-hole symmetry is restored only 
when both antibonding and bonding bands are consid- 
ered. Thus, the maximum covalent effects arise in the in- 
sulating state, in a hybridization gap between two atomic 
levels. The last term would indicate a strong phonon- 
phonon (anharmonic) contribution to the phase transi- 
tions. 

In Appendix A, this analysis is extended to include 
the doping dependence of the LTO transition. Doping 
reduces the interlayer mismatch in two ways, by substi- 
tuting the larger Sr^+ for La^+, and by hole doping the 
Cu02 planes. Insofar as the holes go predominantly onto 
O's, this amounts to replacing 0^~ by the much smaller 
O^^. Thus, in the doped material the layer mismatch is 
reduced, so the tilting instability is driven to lower tem- 
peratures. Detailed calculations based on this picture 
can provide a semiquantitative description of the LTO 
transition, Fig. 2. The data in Fig. 2 come from ex- 
periments on pure Sr doping (summarized by Takahashi, 
et al. |2^), and on Nd substitution studies ||2^. Once 
again, it is found that the transition can be explained 
equally well by either a purely ionic model, or a mixed 
ionic, electron-phonon model. The three frames of Fig. 
2 correspond to the three models of Fig. 3 of Ref. 
with zero (a), weak (b), or strong (c) electron-phonon 
coupling. As in that reference, all three models provide 
comparable fits to the data. 

In principle, the Nd substitution studies could help 
clarify the situation, since Nd is isovalent with La, and 
Nd substitution should lead to only ionic size effects, with 
no doping of the Cu02 layers. However, even this data 
can be fit by any of the above models. 

What is one to make of this? The principle of Occam's 
razor would say that since the purely ionic model (Fig. 
2a) provides the simplest interpretation of the data, then 
in the absence of any other evidence for electron-phonon 
coupling, this ionic model is preferred. On the other 
hand, the cuprates are high-Tc superconductors, with 
considerable evidence for phonon anomalies in the vicin- 
ity of Tc, and with a strong isotope effect when the mate- 
rials are doped away from optimum Tc- In this case. Fig. 
2c can be taken as evidence that the structural transitions 
are compatible with a strong electron-phonon coupling. 
In fact, much more can be said. As the next Section 
will explain, a careful analysis of the purely ionic model 
strongly suggests that there should be strong electron- 
phonon coupling. 
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III. CU-O COMPRESSIBILITY AND STRONG 
ELECTRON-PHONON COUPLING 

A. Doping Dependence of the Lattice Mismatch 

In this Section, the ionic model is analyzed in detail, to 
better understand the doping dependence of the tolerance 
factor. We can begin by relating the interlayer strain r+, 
Eq. 2, to the lattice mismatch. (The middle term of Eq. 
2 has a weak doping dependence, since Uh only varies 
from 1 to 1.2; it can be incorporated into a renormalized 
a)o-) If the equilibrium cell size of the Cu02 plane is 
ai and of the LaO plane 02, then the two layers will 
accomodate to one another at an intermediate cell size 
a, with a2 < a < fli. This will involve a compressional 
strain ei ~ (a — ai)/ai on the Cu02 plane, and a tension 
62 — (a — 0-2) / 0.2 on the LaO plane. The equillibrium 
lattice size a is found by minimizing the elastic energy 



a = 



Ciai + 2C2a2 

Ci + 2C2 ■ 



(3a) 



(36) 



(The 2 in front of C2 comes from the fact that there are 
two LaO layers.) This ensures that the stresses in the 
various layers are balanced, 

Ti - -2T2, 

where Ti = CiCi. Because of the stress-balance equation, 
the only free parameter is T2, or equivalently 02 — a. 
Hence, T^to is a function only of 02 — a. This is conve- 
nient, because a is the measured radius, and the doping 
dependence of 02 can be estimated from tables of ionic 
radii. In contrast, evaluation of ai requires knowledge 
of how hole doping modifies the Cu and O radii. In Ap- 
pendix B, the parameter — 02 is estimated as a function 
of doping, and this result is used to derive the theoretical 
curves in Fig. 2 (see Eq. A2). As a result of the calcula- 
tion, the largest contribution to the doping dependence of 
the transition is found to be the large ionic size change 
associated with 0^~ ^0~, with estimates for the O^ 
radius that are about 14-21% smaller than that of O'^^. 

The same calculation can also explain the Nd sub- 
stitution results. The Nd ion is smaller than the La 
(1.1631 for Nd3+ vs 1.216A for La3+), so Nd substitu- 
tion will worsen the mismatch, and hence should enhance 
the structural transition temperature. This is indeed ob- 
served: there is a strong variation of Tlto with Nd con- 
tent (Fig. 2), which is quantitatively well described by 
the present model. However, the in-plane area of the or- 
thorhombic cell, ab, is virtually unchanged by Nd sub- 
stitution. This suggests that the cell size is fixed by 
the Cu02 planes, while changing the La-Nd size leads 
predominantly to changes in incommensurability strains. 



Correspondingly, it is found that increases in Sr doping 
lead to large decreases in ab. 

This can be understood from Eq. 3b. If Ci >> C2, 
then a ~ fli , so the cell area will be independent of the Nd 
substitutions in the LaO layer, as observed. This is con- 
sistent with the expectation that the LaO layer should be 
more compressible than the Cu02 layer Moreover, 
the Cu02 layer is already compressed via the distortion 
of the CuOg octahedra, with the four planar O's moving 
closer to the Cu, the two apical O's moving away. (This 
is a molecular JT effect, which splits the degeneracy of 
the dri.2_y2 and orbitals on the Cu.) Hence, it will 
be hard to produce any additional Cu02 compression (in 
the absence of hole doping) . 

In order to explain the pressure dependence of the 
LTO transition temperature, dTLxo/dp < 0, Zhou, et 
al. suggested that the Cu-0 bond length is more 
compressible than the La-0 bond length. This result is 
surprising, and is not consistent with the above finding, 
Ci >> C2. Schilling and Klotz |2^] have suggested that 
the anomalous Tlto{p) may be due to pressure induced 
charge transfer. While the sign of the charge transfer is 
the same as that found in other cuprates, direct Hall ef- 
fect measurements have not found evidence for a pressure 
dependence of the hole concentration in LSCO. 



B. Large Electron-phonon Coupling to In-Plane 
Phonons 

The above calculations reveal that the purely ionic 
model is not self-consistent. The dominant factor in 
explaining the Sr doping dependence of the LTO tran- 
sition is the large dependence of the O radius on va- 
lence. However, it is well known in rare earth compounds, 
that such large ionic contractions lead to large electron- 
phonon coupling, accompanied by lattice softening and 
elastic anomalies fl^ . 

When electron-phonon coupling arises through the de- 
pendence of ionic size on valency, a large number of 
phonon modes can couple to the distortion. This occurs 
in other mixed- valence systems as well (see, e.g., Mook, 
et al. [Q). In the present case. Figure 3 illustrates some 
of the Cu02 plane phonon modes which should have a 
strong coupling to the valence fluctuation. Thus, if the 
0~ is localized on a single site, the local contraction will 
involve a Cu-Cu bond stretching mode. Fig. 3a. If it 
can hop between adjacent O's across an intervening Cu, 
it will couple to the 0-0 bond stretching mode. Fig. 3b. 
These modes are expected to be longitudinal within a 
single Cu-O-Cu row on a given plane, but in successive 
rows, the stretch distortions can be either in-phase, as in 
Fig. 3a, b, or out of phase, as in Fig. 3c, d. It turns out 
that, due to vHs coupling, the gaps are larger for the out 
of phase modes. Fig. 3c, d. 

For completeness, certain other modes are illustrated 
in Fig. 3e,f. Figure 3e shows a related, ferroelectric 
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mode, a 2D version of the mode introduced by Egami, 
et al. |]l0|. In addition to the optical modes, there can 
also be coupling to elastic strains (to account for the 
average, homogeneous effects of volume reductions) and 
acoustic phonons (for the local, inhomogeneous effects). 
Fig. 3f. These shear strains function as secondary order 
parameters for the LTO transition. 

Since these are all in-plane modes, it is not clear why 
the actual observed soft mode involves out-of-plane tilt- 
ing. A possible explanation is developed in Fig. 3g. The 
extra hole is assumed to be partially delocalized in the 
x-direction, giving rise to a local shear, shrinking the lat- 
tice a;-axis and increasing the y-axis. Since the layer is 
already under compression, the octahedra rotate about 
the X-axis, to tip the long-bond out of the a — b plane. 
This would explain the highly anomalous feature of the 
LTT phase, that the Cu-Cu distance is shorter along the 
untilted bond than along the tilted one. An analogous 
effect occurs in the LTO phase. 

Recently, Zhou, et al. |^ have also come to the con- 
clusion that in-plane phonon modes play a large role in 
the dynamic JT effect in LSCO. 

Thus, in analyzing the soft-mode instabilities, the fol- 
lowing hypothesis will be adopted. The primary electron- 
phonon interaction involves a planar distortion, but as 
this softens, it interacts with the lower lying LTT-type 
tilt mode, leading to an actual soft mode of mixed sym- 
metry. For a first analysis, this additional complica- 
tion of phonon mode-mode coupling will be ignored, and 
only the planar mode will be included in the analysis of 
electron-phonon coupling. 

C. Other Problems with the Ionic Model 



these two bond lengths are equal. 

This leads to a second problem: on each side of the 
phase transition, the stable phase is the one with the 
longer Cu-0 bond length] This is opposite to the predic- 
tion of the ionic model, which says that the driving force 
for the transition is to reduce the Cu-0 bond length. 

There is also another potential problem. The dot- 
dashed line in Fig. 4 shows the doping dependence of 
the equilibrium Cu-0 bond length, which agrees with 
the measured value in the LTO phase, but changes much 
more rapidly with doping than the experimental value 
in the HTT phase. The break in slope cannot be ex- 
plained by the tolerance factor passing through the value 
unity. Theoretically, this is not expected to happen until 
X ~ 0.42. More importantly, the break in slope occurs at 
different x values at different temperatures. However, the 
HTT data all come from the overdoped regime, x > 0.15. 
It is known ]33| that the uniform phase is metastable 
against phase separation in this regime, and moreover, 
that the doping holes go into additional orbitals - in par- 
ticular, the Cu orbitals. Hence, the change in slope 
may be a signature of this electronic effect - experiments 
at higher temperatures and lower doping are clearly de- 
sirable. Again, this shows that electronic effects must 
be included in order to understand the structural phase 
transitions. 

The dynamic JT model can provide at least a partial 
explanation of these facts. It is important to recall that 
the onset of the tilting instability actually occurs at T*. 
Hence, since both the LTO and HTT phases are dynamic 
JT phases, the HTT-LTO crossover will be controlled by 
long-range strain forces, and need not involve a substan- 
tial change in the Cu-0 bond length. 



There are a number of other features of this structural 
transition which are difficult to understand in a simple 
ionic picture. First, if the driving force for the transition 
is a layer mismatch, with the LaO lattice constant being 
too small, it is somewhat surprising that the La O p lane 
is buckled rather than flat. Billenge and Egami js^] also 
questioned the ionic model, noting that the Cu02 planes 
in Nd2_a;Cea:Cu04 showed local buckling, even though in 
the ionic model, they are supposed to be under tension. 

More importantly. Fig. 4 shows the doping depen- 
dence of the Cu-0 bond length found by Radaelli, et al 
psf . The detailed form of this dependence is inconsistent 
with a purely ionic model. There are actually a number 
of problems. One problem lies in the temperature depen- 
dence at fixed x. If the tolerance factor were the driving 
force for the transition, it would be expected that the Cu- 
O bond length would act like an order parameter, having 
a strong temperature dependence in the LTO phase. In- 
stead, this bond length has a strong x dependence, but 
is nearly independent of T in the LTO phase, whereas 
it has a strong T dependence and weak x dependence in 
the HTT phase. The LTO-HTT transition occurs when 



IV. LINEAR ELECTRON-PHONON COUPLING 
AND THE O-O BOND STRETCHING MODE 

Figure 3 illustrates some of the phonon modes which 
could couple to the large 0-atom valence fluctuations. To 
better understand the nature of the resulting electron- 
phonon coupling, the modification of the electronic band 
structure can be calculated, in the presence of a static 
distortion of the appropriate symmetry. It is assumed 
that the distorted bond lengths lead to a modulation of 
the Cu-0 hopping patameter tcuO, tcuO tcuO i St. 
Assuming the standard three-band model for the Cu02 
planes. Figure 5 illustrates the resulting dispersions, 
showing that several modes can open a gap at the Fermi 
level. The calculations are given in Appendix C. (The 
figure is calculated assuming bare band parameters, t — 
leV, A = 4eV; for the renormalized bands ||3^, the en- 
ergy scale should be reduced by a factor of 4.) 

From Fig. 5, it can be seen that the largest gap is 
associated with the 0-0 bond stretching mode. Fig. 5d. 
The Cu-Cu bond stretching mode has a comparable gap, 
but with a strong fc-dependence, and vanishes at the vHs 
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(5c). This difference is discussed in Appendix C. No gap 
is opened for the FE mode (5e) or the strains (5f), al- 
though in the latter case, the anisotropy splits the vHs 
degeneracy along the X and Y axes. The differences be- 
tween the qi nesting (Fig. 5a, b) and go nesting, (Fig. 
5c, d) lie mainly in the different zone folding associated 
with the doubling of the unit cell. For qi nesting, the X 
point is mapped into F, while M is mapped into Y] on 
the other hand, q^ maps M into F and X into Y . For 
ease in comparison. Figs. 5e,f are illustrated including a 
qg-type zone folding. 

The 0-0 bond stretching mode looks like an ordinary 
ID Peierls instability cum CDW, and has the form postu- 
lated for the vHs-JT effect. Tang and Hirsch had an- 
alyzed the CDW transition several years ago, and found 
that it did not couple to the vHs's. However, they used 
a Cu-only model, so their CDW corresponds to the Cu- 
Cu bond stretching mode, which indeed does not split 
the vHs degeneracy. Fig. 5c. From Eq. C9, it can 
be seen that the 0-0 bond stretching mode. Fig. 3d 
and the shear strains. Fig. 3f, involve a linear electron- 
phonon coupling, as opposed to the quadratic coupling to 
the tilt-mode phonons. This linear coupling will lead to 
large polaronic effects, since the more localized the hole 
is, the closer to 0~ is the corresponding oxygen, with 
a much larger local contraction. Such polaronic effects 
have often been postulated to play a role in the high-T 
superconductivity. 

Direct evidence for anomalously strong elect ron-phon- 
on coupling of the 0-0 bond stretching mode has been 
found by neutron scattering measurements, both in 
LSCO and in YBCO, and closely analogous results are 
found in Bai_^K^Bi03 (BKBO). In undoped LSCO, 
this one-dimensional bond stretching mode is found near 
20THz, with a weak softening near the X-point. Upon 
doping, the softening substantially increases |36) . A sim- 
ilar result is found in YBCO 

In BKBO, it has generally been believed that the struc- 
tural anomalies are associated with three dimensional 
breathing modes. However, in the heavily doped materi- 
als there are soft modes associated with BiOg octahedral 
tilts, similar to those in LSCO jl^. Moreover, no anoma- 
lous softening of the three-dimensional breathing modes 
is found. Instead, the one-dimensional Bi — O bond 
stretching mode softens, again in complete analogy to the 
cuprates . It has proven difficult to model this soften- 
ing in a lattice dynamical model, assuming the local sym- 
metry is cubic. These results suggest that in BKBO the 
structural and superconducting anomalies may be associ- 
ated with a 3D vHs, as in the Bilbro-McMillan model , 
which has been applied to the A15's and BaPbi_a;Bij;03 



V. POLARONIC HAMILTONIAN 

Let us briefly summarize the above discussion. The 
analysis of the doping dependence of Tlto shows that, 
while inter layer mismatch is a dominant driving force, a 
purely ionic model cannot provide an adequate descrip- 
tion of the data. Moreover, even within the ionic model, 
most of the accomodation of the mismatch with doping 
is found to be due to a large variation of the in-plane 
Cu-0 bond length with doping (chemically, with the va- 
lence change 0^~ — >0~). This bond length change itself 
produces a large electron-phonon coupling, with a num- 
ber of in-plane modes. The strength of coupling can be 
estimated from the size of the gap resulting from a CDW 
distortion: the largest gap is associated with the 0-0 
bond stretching mode. Fig. 5d, due to the splitting of 
the vHs degeneracy. However, the symmetry of this dis- 
tortion is not that of the LTO phase, so this distortion 
can only be present in the form of a local, short-range 
order, possibly dynamic. Such short range order has al- 
ready been postulated in terms of the tilting modes 
This picture is supported by experiments, which find sig- 
nificant softening of the relevant modes, and LDA frozen 
phonon calculations, which find a significant phonon soft- 
ening for this mode. 

These results suggest an important modification to the 
earlier calculations of the LTT-LTO-HTT transitions as 
a series of dynamic vHs-JT effects P-|6|,p4[. In these 
papers, the soft mode was identified with the experi- 
mentally observed mode, associated with octahedral tilt 
(plus an associated shear strain). This led to an unusual 
electron-phonon coupling, quadratic in the tilt angle 6, 
and hence in the phonon creation operator - in violation 
of Migdal's theorem. Nevertheless, the resulting coupling 
is estimated to be quite large . The present results 

suggest that the tilting is a secondary effect, and the 
primary driving force is associated with in-plane modes 
(particularly the 0-0 bond stretching mode) that have 
a linear electron-phonon coupling. Hence, there are im- 
portant polaronic contributions, which will be estimated 
below. Note that many of the predictions of the model 
are not sensitive to the particular phonon mode; these 
include the occurence of a pseudogap Q, and of a dy- 
namic JT effect J|6| , which has apparently been observed 
experimentally 0,^. 

In this section it will be assumed that the principal 
electron-phonon coupling is with the planar bond stretch- 
ing mode, and that the tilt mode softens first due to 
phonon mode-mode coupling with these planar modes. 
The complications due to this mode-mode coupling will 
here be ignored, leaving a problem with a linear electron- 
phonon coupling, which can be analyzed using standard 
techniques. It should be noted that the resulting model 
is essentially a vHs-driven CDW instability, and hence 
very similar to a number of earlier calculations |^,^,^ . 
The differences are in the identification of the dominant 
phonon mode (Fig. 5d), and the realization that the 
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physics is dominated by dynamic, short-range fluctua- 
tions, having a different symmetry from the macroscopic 
average. 

A. Starting Hamiltonian 

The three-band Hamiltonian can be written: 

H = He+ Hph + He-ph, (4) 

where the electronic Hamiltonian is 

He - Y.i''i'^% + E 'PP]+SP1+S + E ^^4p,+S + (C-C-)] 
IS s 

+ ^too[pl^^Pj+S' + (c-c.)] + Unji^nji), (5) 

S' 

where {p^) is a creation operator for holes on Cu (O), 
j is summed over lattice sites, S over nearest neighbors, 
S' over next-nearest (0-0) neighbors, and c.c. stands 
for complex conjugate. A = — ep is the Cu-0 split- 
ting, t (too) is the Cu-0 (0-0) hopping energy, and U 
is the on-site Coulomb repulsion. In a slave boson cal- 
culation, the Coulomb term does not appear explicitly, 
but via correlation effects, which renormalize the one- 
electron band parameters. In this case, the eigenenergy 
of the one-electron Hamiltonian, becomes 

Ej: = j + \lijr+^t^sl + sl), (6) 

where A = e^j — Cp, and I have introduced the notation 
(Appendix C) Si = sin{kia/2), Ci = cos{kia/2). 
The pure phonon term is 

+ ]^Y.^l{\,q)Q{\,q)Q{\~q) (7) 

where Q{\,q) {P(X,q)) is the phonon position (momen- 
tum) operator. The form of H^-ph is discussed below. 
In this paper, only a highly simplified model will be ana- 
lyzed, to illustrate the nature of the polaronic coupling. 
Among the simplifications assumed are (1) too =0; (2) 
correlation effects are small enough that they can be in- 
corporated by the use of renormalized band parameters; 
(3) the only relevant phonon mode is the ID 0-0 bond 
stretching mode. 

Technically, assumptions (1) and (2) are incompatible, 
but this is unimportant for the present, illustrative pur- 
poses. Assumption (1) forces the vHs's to fall at half 
filling, where correlation effects are large, and suppress 
the electron-phonon coupling, particularly near the in- 
sulating phase p^ , ^ . In this case, spin-Peierls cou- 
pling should be included, and the resulting analysis is 



not greatly modified. This is the topic of a future publi- 
cation j4|] . In the doped material, correlation effects are 
relatively small, and the present analysis should apply. 
However, a non-zero too is required to move the vHs 
into this regime, and the extra couplings introduced by 
this term greatly complicate the analysis, without signif- 
icantly affecting the underlying physics. 

Assumption (3) restricts the analysis to a single branch 
of phonons, which includes the 0-0 bond stretching 
mode. The distinction between the phonons of the form 
of Fig. 4b and 4d lies in the direction of the phonon 
propagation vector, corresponding to two different cuts 
along this branch. For the cut of Fig. 4b, q = (qi, 
with qi — (7r/a,0) and C £ [0,1]; for cut 4d, q — Cqo, 
Qo = (7r/a,7r/a). Only the latter cut will be analyzed in 
detail. The corresponding phonon operators will be de- 
noted so = Q{^i<l)-< Pq — Pi^^<l)i LOq- = UJo{X,q). 

It should be noted that this same phonon, at if = 0, has 
the character of a FE mode. In a more detailed treat- 
ment of the phonons, this should correspond to a nearly 
pure 0-0 vibration, with the xO's and yO's vibrating 
out of phase. 

The analysis of Section III (Figs. 4,5) suggests that 
the electron-phonon coupling should be included via the 
dependence of t on the ionic separation 

tij ^ t+ (uj - Ui) ■ Vti-j, (8a) 

where the Mi's are the displacements of the ith atoms. 
The Ui can be expanded in terms of the phonon normal 
mode coordinates as 

-^ = E(5^^"-^'Q(^'^1' (8^) 

with e(A, q) the polarization vector and m the ionic mass. 
In the present, single branch approximation, defining 
— —e{q, A) • Vi/ V mN , the electron-phonon coupling 
with the 0-0 bond stretching mode becomes 

He-pH = QMk) (4p,,,r_,- + p1 -rffe), (9) 

with 

9q{k) = 2g^^. (10) 

(To describe the Cu-Cu hopping, it is only necessary to 
make the substitution — > —cosik^ — qx)a/2 in Eq. 10.) 

The resulting Hamiltonian, Eqs. 4-7,9, is almost ex- 
actly of the form analyzed in Ref. |Q (after Balseiro and 
Falicov 1^^) to study the competition between Peierls in- 
stability and superconductivity. The differences are (1) 
the Fermi surface is square at half filling, since the O- 
O hopping term has been neglected (this was done only 
for simplicity, and an appropriate too may readily be 
reintroduced); and (2) the electron-phonon coupling g{k) 
now has an important wave number dependence. 
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B. Link Operator Formalism 

The electron-phonon coupling term can be eliminated 
by a canonical transformation [ [iq|4q -|50|| : 

H = e-'^'^'He'^/^ +^[H,S]. (11) 

In the usual small polaron theory (e.g., Ref. ||5^), this 
transformation must be carried out to all orders in S to 
correctly describe the polaronic band narrowing. How- 
ever, in the present problem, the function S is quite com- 
plicated, due to the 2D nature of the problem. Hence, 
in what follows, I introduce the following approximation 
scheme. I first find an S which exactly eliminates the 
electron-phonon coupling in first order, Eq. 11. I then 
show that in the transformed Hamiltonian iJ, the hop- 
ping parameter is renormalized to i — *■ t{l — S). This 
is precisely the form of the lowest order band narrowing 
found in conventional small polaron theory, and S can be 
shown to have exactly the same interpretation. Hence, 
it will be assumed that the role of higher order terms is 
essentially the same in the present model, to transform 
the renormalization to [\ — S) ^ . The results are 
given in Section V.C.4. 

S can be expressed in terms of the (generalized) link 
operators as 

S^^J2iQ,^AJ:.-^P,-BJ:.), (12) 

where 

k.q / J •' ik.q" i~.k,q^ ^ ' 

2^1,9 



PB±rk,q = plkPyk^q T p^,^_/.fc, im 

These link operators are the generalization of the opera- 
tors introduced by Halperin and Rice into the prob- 
lem of density wave instabilities of an electron gas, and 
first applied to the cuprates by Schulz They also 

bear a close resemblance to the valence bond operators 
studied by Affleck and Marsden . 

A few technical points about these operators should be 
noted. First, the Hamiltonian, Eq. 4, can be rewritten 
in terms of these operators: 



+*.kPl-Ao + %kP3~.k.o)^ (15") 

He-ph = Q?9qi^)Pl+M.r (-^5^) 

k,q 

with t.j: = 2its.j:, s.j: — sin{kia/2). Secondly, the sign 
change in Eqs. 14e-i arises because the hopping term 
in Hf. is imaginary {isx)- Finally, the set of p's is ac- 
tually overcomplete, since, e.g., P2±,k,q = ^Pi±X-q-q- 
Because of the 2D nature of the problem, finding the 
/i's and /li's requires inverting an 18 x 18 matrix. The 
solution is given in Appendix D. 



E ^k,qP^+,k,?^ (13&) C- Polaronic Effects 

1. Effective Interaction 



1=1,9 



Pi±,k,q~ jiPxk-q P'xk-q ^ ^ Oucc the linear electron-phonou coupliug is eliminated, 

we are left with a complicated Hamiltonian involving 
^ ^ electron-electron and quadratic electron-phonon intcrac- 

P2±X,q~ Pj-k'^k-q-^ '^k-qP^i^^ (14^) tious, describing a rich variety of phenomena, including 

phonon softening, structural instability (CDW forma- 
tion), and superconductivity. In this paper, no attempt 

p„, r - = dtp r ^± n^- jIt:, (14c) is made to sort out all of these effects, but only to sum- 

^ yK — q yk — q ^ i-i 1-1^ 

marize the various phenomena which are expected. For 
this purpose, two simpler, approximate solutions are also 
P4-1- J - = P^-^fc -± dt p j:^ (14*^) presented in Appendix D: one is a three-band, but quasi- 

4±,fe,g yk k q k-q yk solution, in which coupling terms involving only the x 

direction are retained. The second approximate solution 

jt J ^ jt J /I /t ^ is an effective one-band solution, in which the three-band 

p. , .- - = rflrfr -=F"- -"r, (14e) , , , 

oiA,« k 1^ q k-q « electronic Hamiltonian is diagonalized, and only transi- 

tions within the antibonding band are retained. The re- 
.j. I suiting equations are qualitatively similar to standard re- 

Pe±,k,q ~ PxkP'-i:k-q~^ Pxk-qPxk^ (14/) sults for polarouic calculatious. However, there are some 
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subtleties, involving band mixing and a potential k de- 
pendence of the couplings. 

The qualitative nature of the resulting polaronic Ham- 
iltonian can be seen by ignoring all the link operators 
except pi in Eqs. 12-14, and approximating 



/ lk,q ' 



A2 + 16*2 • 



(16a) 



(166) 



where AE = Ej:- Ej:^^, ry = 1/2(1 + {A/Atf), and the 
electron-phonon coupling constant can be roughly esti- 
mated as \gq\ Vt/VmN = l3t/{a/2)VrnN, where /? 

is a dimensionless parameter, ~ 3 — 3.5 |^,^. This es- 
timate relies on a simple tight-binding model. A better 
way to estimate g would be via an LDA frozen phonon 
calculation jl^ . 

Using numerical estimates of t = 0.25eV, (3 = 3.5, and 
TiLj^ ~ 50mey, a modest value Ep ~ 6meV is found. 
However, since AE << fito^^ Grilli and Castellani p^ ] 
showed that correlation effects should not reduce g. This 
would correspond to using a bare value t ~ l.ieV in the 
above estimate, leading to Ep ~ 170meV. 



h - ~ -JM 



(16c) 



(see Eq. Dll). The commutator of S with He^ph gives 
the residual interaction terms. 



3. Phonon Softening 



In a mean field decoupling, 
mode softening, 



(2) 



^ leads to a phonon 



— [He-phjS] — He-e + ^.^^^ 



e—ph ' 



with (Eq. D9) 

He-e = -Xl''lfe^«-^'^(fc)Pl+,A:,g-Pl+,fe^g- 



(17) 



(18) 



k,k' 



(2) 

The remaining term, H^J^^f^ is complicated in the three- 
band model, because it mixes in different link operators, 
p (see Eq. D13). However, it simplifies in the one-band 
model, to 



H. 



(2) 
—ph 



J2 ^^^^[«fc+g-4+2/fc " °'k-q^l+g'f'k-q]^ In Eq. 22, the E's can be replaced hy E' ^ E - Ef 



2 
Wo 



nr: ^MEr-Er -.) 
k-qo'^ k k — qo' 



Er -)^ 



(21) 



where —< ip'^-ip^ >, and the precise form of 7^ is not 
important for present purposes (see Eq. D22). In a self- 
consistent calculation, the bare wq in the denominator of 
Eq. 21 would be replaced by luq. Making this substitu- 
tion, a phase transition occurs when lvq is renormalized 
to zero, or 



7fc 



k-qo 



Ek- 



Er - 

k-qo 



(22) 



(19) 

where ij: is the wave function of the antibonding band, 
/i (Eq. D20) is closely related to f^-j: above, and is 
given by Eq. D18b. 



Now if the Fermi energy falls at the vHs, then EL ^ 
— i?^, so the transition temperature Tq is given by 



^0 =< Ik > 



-B 



N*{E)dE 
2E 



tanh{ 



E 



2kBTa ' 



(23a) 



2. Polaron Binding Energy 

The polaron binding energy can be defined from Eq. 
18 as 

Ep = liniAE^Q fei k.-q9q{^) 



\^ \9q{k)\'^V 



2+2 



(20) 



where < 7^ > is the average of 7^: over angles, B is 
the half bandwidth, and N*{E) is an effective density of 
states (dos) 



iv*(i?) = -E 



7fc 



<lk> \QEldk.,\ 



(236) 



The integral in Eq. 23 must be evaluated numerically, 
but certain general features can be directly extracted. 
The result is sensitive to the angle dependence of 7^. If 
7^: is a constant, then the results are the same as found 
earlier N*{E) is the actual dos, which diverges log- 
arithmically at the vHs, greatly enhancing To- This ap- 
pears to be consistent with the gap found for the 0-0 
bond stretching mode. Fig. 5d. However, the approxi- 
mation of retaining only one pair of link operators is too 



9 



crude to reproduce this feature, and the model 7's are 
found to have a more complicated angular dependence. 
Indeed, it is expected, on analogy with the results of Ap- 
pendix C, that extreme care must be taken to derive the 
correct angle-dependence of the gap. 

For other phonon modes, however, is expected to 
have an important angle dependence. Thus, for the 
Cu-Cu bond stretching mode, 7^ ~ jQsin^kxa (see the 
discussion below Eq. C8). In this case, N*{E) ~ 
16-\/-B^ — E'^ /{ttBY, remaining finite at the vHs, E — G, 
and 



B 

\.\A—e 
2 



l/A 



with 



(24a) 



(246) 



with TVi = 16/7r^i?. Thus, for the Cu-Cu bond stretching 
mode, the CDW has a logarithmic instability, and not the 
It? instability expected for a vHs-related singularity, due 
to the fact that the CDW gap for this mode vanishes at 
the vHs, Fig. 5c. This is consistent with earlier work 



4- Polaronic Band Narrowing 
(2) . 

The term H^J_, also induces a polaronic reduction of 
the bandwidth Bfl]. In the Appendix, this reduction is 
calculated to lowest order in 5 as i ^ t{l — S), Eq. D15, 
with 



1. 



where the Bose function is 



-^9 ghuig/kBT _ I ■ 



(25a) 



(256) 



In all of the above expressions, the Fermi energy enters 
only through the Fermi-Thomas function, nj:. The pola- 
ronic band narrowing, Eq. 25, however, is independent 
of uj:, and hence of Ep\ It depends only on the details 
of the band structure; in particular, there is a significant 
band narrowing only when 



(25d) 



When this condition is satisfied, the bandwidth reduction 
is similar to the usual phonon renormalization, 1/(1-|-A), 
with A ~ e'^. 

The condition, Eq. 25d, is discussed in more detail 
in the following section. It is found that this condition 
can be satisfied only in the immediate vicinity of the vHsl 
This may explain the presence of 'extended vHs's' |Q in 
the cuprates. Figure 6 illustrates a calculation of the 
modified band structure of Bi-2212, if it is assumed that 
there is a phonon renormalization of the energy bands 
by a factor A = 3 within u}q = AOmeV of the Fermi 
level. The value A = 3 is in satisfactory agreement with 
the renormalization factor exp{S) found above. More- 
over, a similar band narrowing is found near the vHs in 
YBa2Cu408 0, even though the vHs is ~ 19meV away 
from the Fermi level. 

If terms in too are retained, the parameters t and too 
are renormalized by different factors, so polaronic effects 
could lead to a change in the shape of the Fermi surface. 
Analogous effects found in the f-electron problem js^ j 
suggest that direct 0-0 hopping may be more seriously 
affected than Cu-0 hopping. 

In ordinary small polaron theory | ]50| , the higher order 
commutators play two roles: for the coherent hopping, 
they lead to a polaronic band which becomes increasingly 
narrow as temperature is raised, as in Eq. 25. However, 
they also lead to an incoherent hopping due to nondiago- 
nal transitions (in which the number of phonons bound to 
the electron changes during the hopping process). This 
latter effect is not accounted for in the present calcula- 
tion, which treats S to lowest order. 



In conventional polaron theory, inclusion of higher order 
commutators in S leads to an exponential bandwidth re- 
duction. Hence, I hypothesize that the main effect of 
including higher order commutators with S in Eq. 11 
would be the renormalization 

t te~^. 

However, since S can have either sign, a more appropriate 
approximation would be 



5. CDW vs. Superconductivity 

The CDW and superconducting instabilities can also 
be studied in terms of He-e- For the CDW instability, 
only the term in q = qo is retained |^,^. With this 
approximation, i?e-e becomes, in mean field. 



CDW 



= -Go ^ Cx 



l+,k,qo ■ 



(26) 



t 



2t 



1 + • 



(25c) 



By comparing Eq. 25a and Eq. 20, it can be seen that at 
low temperatures {Ng^ ~ 0) 5 ~ Ep/2hijj^, so using the 
above estimate for Ep, exp{S) ~ 5. 



Since _ffe + Hcdw is quadratic, it can be diagonalized 
by the Bogoliubov-Valatin transformation. The resulting 
dispersion satisfies Eq. C6, if Gq = 25t. Again, the 
present approximate models lead to a Go with a more 
complicated angular dependence. 
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For superconductivity, those terms in Eq. 18 are re- 
tained, for which k' = —k, and a' = — ct, where the 
cr's are the spin indices, which had been neglected up to 
now. This leads to a conventional, vHs-enhanced super- 
conductivity 1^, except that the gap equation contains 
an angular factor similar to 7^! Such an angular factor 
can, in principle, lead to a non-s-wave gap. 

VI. DISCUSSION 

A. lonic-Covalent Crossover and the vHs 

In their discussion of electron-phonon coupling induced 
by valence fluctuations in f-electron metals, Sherring- 
ton and Riseborough (SR) introduced a Hamiltonian 
with a somewhat different form of He-ph- SR assumed 
a direct coupling to the valence fluctuation, which in the 
present context would give 

^f-1..-E5.rP^,.P.>Qo' (27) 

Here, while the tight-binding electron-phonon coupling 
has a different form, Eq. 9, the resulting Hamiltonian 
automatically includes a term like Eq. 27, but with 
quadratic phonon coupling, Eq. D13. 

The present calculation sheds some light on a feature 
that SR were attempting to describe: the competition 
between covalency and ionicity. They pointed out that 
earlier theories of electron-phonon coupling in f metals 
ranged from strongly covalent |]57| , with each rare earth 
ion having average values of ionic size and f-electron oc- 
cupation number, to strongly ionic [Q, with well-defined 
/" and /"~^ ions. SR suggested that this competition is 
represented by terms describing hopping, t, and electron- 
phonon coupling, He-ph'- He-ph << t represents the co- 
valent limit, where the valence can freely fiuctuate, while 
the opposite limit, H^-ph >> ^7 is the ionic limit, in 
which each ion has a well defined valence. 

This can be restated: the strong electron-phonon cou- 
pling is associated with valence fluctuations: a hole hops 
onto an 0^~, causing its radius to collapse, leading to 
an enhanced local value of t. For the collapse to occur, 
the phonon must be able to follow the electronic motion, 
<^ph ~ {hopping rate). For an 'isolated CuOe molecule', 
this hopping rate is ^ t/h >> ajph- Near half filling, 
strong correlations reduce the hopping, but the inequal- 
ity still holds, J ~ ISOmeV > hojph ~ 30 — 70meV. 

In a band, however, the electronic motion is restricted 
by the motion of other electrons, and the relevant fre- 
quency becomes q- vf, where q is the phonon wave num- 
ber and Vf the Fermi velocity. At an arbitrary point on 
the Fermi surface, this quantity is of comparable mag- 
nitude to t/h, but Vf vanishes at a vHs\ Figures 7a, b 
shows a map of the vf surface for a variety of different 
Fermi levels, while Fig. 7c (solid line) shows the mini- 
mum value oiqQ-VF, as a function of Fermi energy. What 



is plotted is the velocity normalized to a zone boundary 
phonon, v* = hirvFi/a; for a phonon away from the zone 
boundary, the curves should be scaled by the ratio qa/ir. 
(In this plot, the renormalized band dispersion parame- 
ters are used, t = 0.25eV, A — leV.) From these figures, 
it is clear that there is only a narrow range of energies for 
which q-VF^ t^ph , particularly near the zone boundary. 
Moreover, the slowest velocities are associated with hop- 
ping in directions near X and Y - i.e., towards the Cu's. 
For these directions, and sufficiently close to the vHs, 
the electronic motion is slow enough that the molecular 
distortions can keep up. 
The condition 

^ph — q-VF (28a) 

can be reinterpreted in a number of ways. First, for 
acoustic phonons the sound velocity is Cg — ujph/q, and 
Eq. 28a becomes 

Cs^VF (285) 

- the 'surf riding' condition that signals the onset of Lan- 
dau damping. Finally, for large g, the electronic disper- 
sion should be treated more carefully, and Eq. 28a should 
be replaced with 

for some electronic wave number fc, where E-j^ is the elec- 
tronic energy dispersion. Equation 28c should be com- 
pared to the resonant denominators of Eqs. 16a and 21. 
The right-hand side of Eq. 28c is plotted as the dashed 
line in Fig. 7c, and it can be seen that use of this con- 
dition broadens the range of energies over which the in- 
equality is satisfied. This last equation is the condition 
for a resonant electron-phonon interaction - i.e., both 
energy and momentum are conserved. Such a coupling 
will lead to large modifications of both the electronic and 
phononic dispersions. 

This suggests a new interpretation for the vHs. Away 
from a vHs, electronic hopping is very fast, and it is 
hard to say just which atom the electron is on instanta- 
neously. The band is a mixture of Cu and O, with both 
at an intermediate valence - i.e., the coupling is cova- 
lent. The Born-Oppenheimer (adiabatic) approximation 
is valid: the electronic motion is so fast that it responds 
instantaneously to any nuclear motions. Near a vHs, the 
opposite, anti-adiabatic limit is valid. The electronic mo- 
tions are slow compared to the nuclear motions - for at 
least some of the electrons, vf = 0. In this limit, the 
valence of an individual atom is much more definite. For 
instance, in the fast hopping regime, the local distortions 
associated with vs do not have a chance to de- 
velop, so all O's have a common, intermediate value of 
radius. Near the vHs, on the other hand, Cu-0 hopping 
is slowed down, allowing local distortions to develop, so 
that an 0~ can be clearly distinguished from an 
by its smaller radius. In short, near a vHs the nuclear 
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motion can 'instantaneously' follow the electronic mo- 
tion, and the bonding goes from covalent to substantially 
more ionic. This leads to the strong enhancement of the 
electron-phonon coupling. 

Since correlation effects also act to slow down the Cu-0 
hopping, this may be an explanation for why correlation 
effects act to pin the vHs nearer to the Fermi level at half 
filling. 

B. Future Work 

The present paper is the second of a projected series of 
papers whose aim is to provide a firm basis for analyzing 
the vHs-JT model within the three-band model of the 
cuprates. The first paper (XA) jsij had two purposes: 
the main purpose was to extend the slave boson calcu- 
lation of correlation effects to the full 3-band model; the 
second, to explore a possible linear electron-phonon cou- 
pling induced by spin-orbit coupling. The present paper 
(XB) further explores electron-phonon coupling in the 
Peierls regime, when correlation effects are expected to 
be weak. A strong, linear component is found to con- 
tribute to the dynamic JT fluctuations, leading to im- 
portant polaronic effects. 

A related paper (XIA) presently under prepa- 

ration, analyzes electron-phonon coupling in the spin- 
Peierls regime, where correlation effects are strong - in 
the 2D Heisenberg limit near half filling. In future pa- 
pers, I hope to extend these results to the three-band 
model. At this point, I should be able to study both 
Peierls and spin-Peierls regimes in parallel. Only then 
will I be able to give a detailed calculation of the dop- 
ing dependence of the structural transitions, including 
the 'pseudogap', and to better understand the possible 
nanoscale phase separation in these materials. 

A brief mention should be made about spin-orbit cou- 
pling, discussed in XA. Spin-orbit coupling leads to inter- 
esting modifications of Fermi surfaces, and, since it leads 
to a linear electron-phonon coupling, it could play a role 
in driving structural instabilities. This effect competes 
with the dynamic JT effect (a competition well known 
in molecular JT problems). It has been suggested that 
spin-orbit coupling is too small to drive any structural 
instabilities. However, a simple estimate (Appendix E) 
shows that, near a vHs, this suggestion may not be cor- 
rect. 



VII. CONCLUSIONS 

The present calculations clarify the nature of electron- 
phonon coupling in the HTT-LTO-LTT phase transitions 
in LSCO, LBCO. First, the dominant effect in the transi- 
tion is the change of Cu-0 bondlength with doping. This 
large valence dependence of the oxygen ion radius is a sig- 
nature of a strong, linear electron-phonon coupling. How- 



ever, as the corresponding phonon modes soften, they in- 
teract with the very soft tilt modes, driving them unsta- 
ble. Thus, although the dominant coupling is linear, the 
final instability involves quadratically coupled phonons. 

Therefore, the full theory involves both holes and two 
types of phonons - the stretch modes and the tilt modes. 
However, the main effects on the holes are produced by 
the stretch modes, and including only these leads to a 
major simplification within the vHs-JT model, in that 
the earlier, quadratic phonon coupling to electrons is re- 
placed by a linear coupling. This linear, in-plane coupling 
can now be treated by conventional techniques, provid- 
ing a possible microscopic basis for the many polaronic 
models of superconductivity |5^] . A striking prediction of 
the model is a significant band narrowing, only in the im- 
mediate vicinity of the vHs, thereby providing a plausible 
explanation for the existence of 'extended vHs'. 

The reason the band narrowing is localized near the 
vHs is connected with a special property of the vHs. Near 
a vHs, the material is in the 'ionic limit', in that the 
electronic motion is slow compared to phonon motion, so 
valence fluctuations are well defined. 
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APPENDIX A: DOPING DEPENDENCE OF LTO 
TRANSITION 

Figure 3 of Ref. jsj showed that the phonon softening 
associated with the LTO transition in La2Cu04 could 
be modelled in terms of a microscopic picture involv- 
ing phonon anharmonicity and electron-phonon coupling. 
The data were insufficient to fix values for all the param- 
eters, and in particular the electron-phonon interaction 
parameter could be varied over a large range, including 
zero. The same model can be applied to explaining the 
doping dependence of the transition temperature, T^tO: 
for either Sr or Nd doping. Figure 2. 

From Eq. 35 of Ref. Tlto is the solution of 

oj^ + [gPi + Pa - a^iV(0)]Ai = 0. (Al) 

Here, wq is the 'bare' phonon frequency, the P^'s rep- 
resent anharmonic (quartic) phonon-phonon couplings, 
tte is an electron-phonon coupling parameter, and Ai cx 
coth(uJi/2kBT) is proportional to the fluctuations of the 
tilt amplitude (the frequency tui is defined in Ref. [^). 
In the absence of electron-phonon coupling, Eq. Al has 
a solution only if ujq is negative. Such negative values 
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can be caused by strains due to interlayer mismatch, and 
to explain the doping dependence of TltOj it will be as- 
sumed that Eq. 2 holds, which can be rewritten as 

uol^ul- f (a - as). {A2) 

(Since rih only varies from 1 to ~ 1.2, the main effect 
of the middle term in Eq. 2 can be absorbed into a 
renormalized value of wq-) 

The temperature dependence of Eq. Al comes from 
that of Al and that of -/V(O), whereas the doping (x) 
dependence is assumed to come from that of the lattice 
mismatch, (a — 02), and that of N{Q). The lattice mis- 
match is estimated in Appendix B. The hardest factor 
to quantify is the dos, iV(0), due to the important role 
of correlation effects. In the absence of correlation ef- 
fects, assuming a rigid-band doping dependence and a 
vHs shifted away from half filling, iV(0) can be written 
as 

jV(0) = --L/n( ^ =), (A3) 

where B is the bandwidth, and /i(x) is the distance in en- 
ergy between the vHs and the Fermi level. This equation 
should hold in the doped materials near optimum Tc, but 
near half filling correlation effects become important. 

Corelation effects reduce the electr on-p honon coupling, 
Oe, for the usual Peierls interaction but there is 

also a spin-Peierls coupling, for electron-phonon coupling 
stronger than a critical threshold |3^,^. The strength 
of this effect is comparable at half filling and at the vHs: 
while nesting is better at half filling (the pseudofermions 
have a square Fermi surface), vHs nesting is suppressed 
by matrix element effects E^]. Hence, the conventional 
nesting at half filling will be of comparable strength to 
the vHs nesting in the doped material. Therefore, for the 
calculations of Tlto in Fig- 2, a simplifying assumption 
is made - that the vHs is pinned to the Fermi level, but 
that the strength of coupling (represelted by N{0)) is 
independent of doping. 

In analyzing the T-dependence of the soft mode fre- 
quency in the undoped La2Cu04 it was found that 
the transition could be equally well fit by either a purely 
anharmonic model or one with a strong component of 
electron-phonon coupling. In analyzing the doping de- 
pendence, it is found that the same ambiguity of pa- 
rameter values remains, and Fig. 2a-c correspond to es- 
sentially the same parameters as Fig. 3a-c of Ref. 
Identical parameters cannot be used, since the transition 
temperature at zero doping is considerably higher for the 
present data than for the earlier (the fits in Ref. |^ found 
Tlto between 455-508K). There is no unique technique 
for adjusting the parameter values to the new T^to with- 
out fitting the entire transition, for which data are not 
available. However, Tlto is directly related to the 'bare' 
phonon frequency, wq, so this parameter was increased 
slightly (by '-^ 4 — 8%) to correct for the higher Tlto val- 
ues (corresponding to fit values of Tlto — 540 — ^7&K 
at X = 0). 



The parameters iwo, wo; and f are listed in Table I. 
The other parameters are the same as in the respective 
parameter sets of Table II of Ref. |^. Note that the 
lattice mismatch is clearly responsible for the instability: 
whereas lo^ is negative, signalling an instability, the bare 
parameter Cjq is large and positive. 

One feature should be noted: the renormalized Nd pa- 
rameter, (a — 02)* from Appendix B, had to be used, to 
get good fits. This is directly related to the fact that 
the electron-phonon coupling term was assumed to be 
doping independent. If a^N{0) is assumed to be smaller 
at X = 0, then electron-phonon coupling would favor a 
higher Tlto in the doped material, a trend that would 
have to be compensated by a larger value of T. To ex- 
plain the Nd-doping results then would require a smaller 
value of (0 — 02) when y is varied. Contrawise, if agiV(O) 
increases as half filling is approached, a larger value of 
the y-coefficient of (0 — 02) is needed. This result can be 
turned around: if we knew perfectly the theoretical form 
of (0 — 02), we could determine from experiment the pre- 
cise doping dependence of the electron-phonon coupling 
coefficient. Unfortunately, the uncertainty in parameters 
is too large to determine whether alN{0) increases or 
decreases with doping. 

APPENDIX B: ATOMIC RADII 

In testing an ionic model for the LTO transition, it is 
necessary to estimate the effective ionic radii in the var- 
ious layers as a function of doping. This will be done in 
two ways: by using tabulated values of radii |]60[| and via 
valence bond sums pi] ]. Exact agreement with experi- 
ment should not be expected, since the tabulated values 
are averages over many compounds. Nevertheless, an ap- 
proximate agreement is found. 

1. On the LaO Plane 

Commensurability strains will arise in LSCO whenever 
the tolerance factor, Eq. 1, is different from unity. It 
should be noted that, in this equation, both dLaO and 
dcu02 are functions of doping. In undoped La2Cu04, 
the Cu02 planes are under compression, t < 1. The 
atoms in the LaO layer are easier to discuss, since they 
can be treated in the ionic limit. The radii will first 
be estimated from Shannon's tables [Q, assuming the 
La site to be 9-fold coordinated |6^. In this case, as- 
suming the conventional O^" radius, rQ2- — 1.40A, the 
relevant radii are rLa^+ = 1.216A, rpf^3+ — 1.163A, 
rBa^+ = I.47I, rsr2+ = I.31A, and rca^+ = I.ISA. 
Thus, for undoped La2Cu04, dLaO — 2.616A, while in 
an Sr and Nd co-doped material, La2-a;-ySr2:Ndj^Cu04, 
dLao = (2.616 + 0M7x - 0.027y)A. Thus, Sr doping re- 
duces the interlayer strains, while Nd substitution makes 
them worse. 
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In the T' structure, the rearrangement of the apical O's 
to interstitial positions expands the lattice, and leaves the 
Cu02 layers under tension |^ , and leading to materials 
which can be electron doped. These structures will not 
be further discussed here. 



2. On the Cu02 Plane 

Shannon's tables list rcu^+ = 0.77A, rc„2+ = 0.73A, 
and rc,(3+ = 0.54A. Assuming the Cu is all in the 2+ 
state, this yields dcuOi = 2. ISA, to be compared with 
02 = dLaol^f^ = 1.850A in La2Cu04. This severe mis- 
match is largely corrected for by a large molecular LT 
effect in the CuOg octahedra: the electronic degeneracy 
of the Cu d levels {dx^^-y^, dz^) is split by a molecu- 
lar distortion in which the planar CuO distance becomes 
much smaller than the apical distance. A rough esti- 
mate of the effect can be made from the Shannon val- 
ues. For Cu^+, the one hole in the Cu d-shell is gen- 
erally an equal mix of the dx2_y2 and the d^^ orbitals. 
In the JT state, the hole is completely taken from the 
dx2_y2 orbital. Hence, for in-plane lengths, the relevant 
orbital has the same occupation as in the isotropic Cu'^+ 
ion, with a corresponding length dtcuO^^"^ — l-^^A. 
For the apical lengths, the dz2 orbital has the same 
length as in Cu^+, or d'o'^Q^ ~ 2.17A. Experimentally, 
these lengths are found to be d'^'Pl''"^ ^ 1.9051 and 
^apical ^ 2.4211 in La2Cu04 at lOK. For the apical O, 
the agreement could be further improved by noting that 
screening is very weak along the c-axis, so the Pauling 
value ||,p§ rc„i+ = 0.961 is more appropriate, yield- 
ing d^™^ 2.361. 

To extend the above picture to the doped regime, it 
is necessary to understand where the doped holes go. 
Experimentally, it is found that hole occupation of the 
Cu dz2 orbital only becomes apparent for x > 0.15 - 
i.e., just above the optimal doping for superconductivity. 
Hence, in the low-doping regime it can be assumed that 
the holes go primarily onto the in-plane O's. Chemically, 
this would correspond to the appearance of the species 
. This is a rare species, and its ionic radius is not 
well established. The present analysis will provide an 
estimate of its value. 



so 



0-02 = (0.051 -0.118a; + 0.019y)l (S3) 

at room temperature. Note that there are additional 
small corrections, which are neglected. In particular, no 
correction is made for the fact that the La and the O are 
not strictly coplanar (and similarly for the Cu and O). 
Also, there is no temperature dependence of the lattice 
constants, and all 9 La-0 distances are assumed to be 
the same. 

The above model predicts too small a change of a — 02 
with y to agree with experiment. Empirically, the ratio 
of the coefficients of y and x can be estimated from the 
change in Tc with these parameters. Now 



Tlto = To - T^x + Tyy, 



m 



with To ~ 546K, ~ 2456i^ (from Ref. l65|), and Tj, ~ 
560K (Ref. H). This yields T^/Ty = 4.39, compared to 
the predicted value, 0.118/0.019 = 6.21. In the following 
subsection, we will see that this lattice mismatch can also 
be estimated from valence bond sums, and here the y- 
coefficient is found to be too large. In the present paper, 
an effective average value will be assumed, to give better 
agreement with experiment: 



(a - 02)* = (0.051 - 0.118a; + 0.027j/)l. 



{B5) 



Note that the calculated transition (Fig. 2) terminates 
at a doping x ~ 0.21, where lUq becomes positive (Eq. 2), 
while the interlayer strain vanishes only when a — a2, at 
a much larger doping, x = 0.43 when y = 0. 



4. Valence Bond Sums 

In view of the importance of the above results, it is 
preferable to recalculate 02 by another technique. The 
method of valence bond sums ||6l|| has been applied by 
several authors to the high-Tc superconductors [|66|-p8[. 
In this model, the valence associated with a bond between 
atoms A and B is written as 



/TABO-rAB)^ 

Vab = exp[ 

ro 



{B6) 



3. Doping Dependence of Tlto 

Using the above estimates, in La2Cu04, 02 — 1.850l, 
oi = dcu02 = 1-941, and a = 1.9051, for a reasonably 
consistent picture. The doping dependence of 02 is 

02 = (1.850 + 0.033a; -0.019y)l, [Bl) 

while, experimentally 

a = (1.901 - 0.085a; + 0.0?/)l {B2) 



where tq = 0.371 is a constant, tabg is tabulated for over 
a hundred possible pairs [|6ll , and tab is the measured 
AB distance. The total valence of atom A is then 



m 



where the sum is over all close neighbor atoms, Bi. This 
can be rearranged to give the ideal radius for a given 
atom: 



TAB = TABO + 'roln[—-), 
VA 
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where Na is the coordination of atom A {N^a = 9). 
Taking VLa = VniI = 3 and Vsr = 2, and using tabulated 
values |6|] for tabo, it is found that 

a2\vai.bond = (1.823 + 0.069x - 0.0472/)!, (B9) 

so 

{a-a2)vai.bond = 1 .31 (0.060 - 0. 1 18x + 0.036y) A. (BIO) 

The coefficient of y is here found to be much larger than 
that estimated from the Shannon tables, and even larger 
than the experimental value. 

5. The Radius of O" 

Now, from Eq. 3b, 

2C 

{ai-a) = —^{a-a2). {BU) 

Since oi should be independent of y, the experimental 
observation that a is approximately independent of y re- 
quires that C2 << Ci. For example, if 2C2/C1 ~ 0.1, 
then a would have to include a term ~ 0.003y. This 
would then yield 

ai ~ (1.904 - 0.097x)A. (S12) 

(Note that, at x = 0, this yields a smaller value of oi 
than the chemical estimate, oi ~ 1.94A, made above.) 
Combining this with Eq. Bl yields a tolerance factor 

i = — ~ 0.972 + 0.069a; - O.OlOy, 

ai 

or t = 1 at a; = 0.42, when y = 0. 

If all the doping is assumed to go into the oxygens, 
then setting x = 2 is equivalent to converting both pla- 
nar O's to O", yielding rp- = (1.40 - 0.194) = 1.2061. 
Grenier |6^ has postulated that O^^ is too large to read- 
ily diffuse, and that the diffusing species in YBCO and 
the other cuprates is probably O". He estimates its ra- 
dius to be ~ 1.1 A, comparable to the present estimate. 
From the valence bond sums, a similar estimate can be 
made, rp^- — tq- ~ roln2 ~ 0.2561. 

The large increase in radius on going to is consis- 
tent with the suggestion that is inherently unstable, 
being stabilized in a lattice by the Madelung potential 
jTOf . This near instability leads to a large, nonlinear po- 
larizability which has been hypothesized to lead to ferro- 
electricity in the perovskites ||70[|. 



APPENDIX C: IN-PLANE MODE COUPLING IN 
A TIGHT BINDING MODEL 

In this Appendix, the electron-phonon coupling as- 
sociated with the various structural distortions will be 



estimated in a tight-binding model, assuming that the 
coupling arises from the bond-length dependence of the 
hopping parameter, tcuO- In the standard three-band 
model, assuming too ^ for simplicity, a very gen- 
eral case can be handled. Assume that the distortions 
have led to a \/2 x unit cell, with alternating hop- 
ping parameters ^i, ^2, ^3, and along the a;-axis (Fig. 
3a, c), and t along y. Then four possible distortion modes 
are (1) Cu-Cu bond stretching mode, ti =^ t2 ~ t ~ 6t, 
^3 = ^4 = t + St; (2) 0-0 bond stretching mode, 
ti = t4 = t — 5t, t2 = tj, = t + St; (3) ferroelectric 
(FE) mode, ti = t^ = t - St, t2 = t^ = t + St; (4) shear 
strain, ti = t2 = t^ — t^ — t — St. For the shear strain 
mode, the y-hopping must be reset to ty = t + St, and for 
comparable distortions, St must be smaller by a factor of 
than for models 1-3. 
For these distortions, the energy dispersion is given by 
the solutions of 

iW-Wj:){W-Wj:_.) = T, (CI) 
with W = E{E~ A), 

y _ t^ 

Wj: = At'' si + £^l^h±± _ + <3t4)c., (C2) 

r - liitj - tlf + {tl - tlf + 2{t^ti - t2t^f 

+2(tii2-t3i4)'(l-2c2)], (C3) 

and Si = sin{kia/2), Ci = cos{kia/2), Cx = cos{kxa). 
The values of are also given by Eq. C2, with the 

substitution k k — qin the arguments of the trig func- 
tions. In particular, for q = qo (Fig. 3c-f) Si — > — Ci, 
Ci ^ Si, i = x,y. The same dispersion relations also hold 
for the distortions of Fig. 3a, b, but with q — qi. In this 
case, the terms in Sy remain unchanged in W^_^. The 
calculated energy dispersions are illustrated in Fig. 5, 
assuming t = leV, A = AeV, and St — O.leT^. Note that 
the Cu-Cu and 0-0 bond stretching modes give identi- 
cal dispersions for a qj distortion (Fig. 5a, b - this can 
be shown analytically from Eq. CI) but very different 
dispersions for the go distortion (Fig. 5c, d). 

In Eq. C3, r is the square of an umklapp scattering 
gap. For the four various modes, it can be simplified to 

f Git'^St'^slcl if i=l: 
r, = <^ m^t^ if i=2; (C4) 

I if i=3,4. 

Note that there is no gap for the FE mode, but a large gap 
^ St opens up for the other modes, although it vanishes 
at the vHs for the Cu-Cu bond stretching mode. When 
7 = 0, the dispersion becomes E — E-^, with 
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For r ^ 0, by defining F{E) = E ~ A/2, Eq. CI can be 
rewritten 



F\E^) + F^{E~^J 



Mi 



FHEj:)-F^E^^_.) 



)2 + r. 



iC6) 



Explicitly, for either the 0-0 bond stretching mode or 
the EE mode, 



for the Cu-Cu bond stretching mode: 

and for the strain: 

Wj: = 4(t - St)'^sl + 4{t + St)'^s 



2 „2 

y 



{C7a) 
iC7b) 
(C7c) 



When Ej^ — -^fc-g- l^-rge, the dispersion is essentially 
given by Eq. C5. However, when Ej: = Ej:_^, the exci- 
tation spectrum is gapped. Eor the 0-0 bond stretching 
mode, 



E^E^,± 



2tSt 



(C8) 



to lowest order in St. Since + Sy = 1 everywhere on 
the Fermi surface, the gap is isotropic, as seen in Fig. 
5d. Similar equations hold for the remaining two distor- 
tions, except that the factor 2t in the numerator becomes 
StSxCx for case 1 (Cu-Cu hopping) and 2-\/2i(s^ — Sy) for 
case 4 (shear strain) . Even though there is no gap in the 
FE band, the distortion still leads to a net energy low- 
ering, but it is ^ St^. This is a 2D version of the model 
discussed by Egami, et al. ||ic[] . 

The gap at the vHs can be calculated as 



AE^ = 



2(3,St 



i — 1,4 labelling the four modes, with Pi — P3 = 0, l32 — 
1, and /34 — \/2. For the strain mode, this represents a 
splitting of the vHs degeneracy at the X and Y points. 



APPENDIX D: POLARONIC CALCULATIONS 

1. Canonical Transformation: Full Three Band 
Model 

The complete, fully 2D solution to Eqs. 11-14 involves 
a system of 18 equations in 18 unknowns, the /'s and /I's. 



Ten unknowns can immediately be eliminated by direct 
substitution. Defining = fi/'-Oq, these are 



ha = 



hn = 



he = 











hujq 




%kh + %k^ 


-g-/4 






Kkh + t.k- 


-q/2 



/15 = - 



\^xk-qf'^ + ^xfc/2 + tyk-^qfs + ^ykfil 

huj^ 



[Dla) 



(Dlb) 



(Die) 



(Did) 



(Die) 



with the five equations formed by interchanging hi <-> fi. 

Solution of the remaining eight equations is more in- 
volved. Defining column matrices 



H 




{D2a) 




{D2b) 




then the remaining equations become 
AH + BF ^ G, 



AF + BH ^ 0, 



with solution 



H^-[{A + B)-' + iA-B)-']G, 



F=-[{A + B)-^ -{A- B)-^]G. 



{D2c) 

(DSa) 
{D3b) 

(DAa) 
{D4b) 
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Introducing the notations a = t^^_^, b = t^^, c = tyj^_^, 

d = ty^, and e 
can be written 



d = ty^, and e = fuj^A, the matrices A = PRP and B 



B 




{Dha) 














/a 














b 














c 





Vo 








d 



{D5b) 



R 



a 


2 


1 


2 


2 


f3 


2 


1 


1 


2 


7 


2 


2 


1 


2 


5 



with 



a = 



a-- 



7 = 



(5 = 



VF-a2 



(£)6a) 

(I?6c) 
(£)6rf) 



and W = n^ujl+o? + b'^ + c^+ d^. 

Defining H± = {A±B)-^G, W± = W±e, A = c?^(?, 
B = b-^ + d"^, the explicit solution is 



■ 



{W+W- - 4AB){W± - A)' 
hi± = W^{W± -A) + c^{W^ - 4B), 
h2± = -2ab{W±-A), 
h3± = -ac{W^-4B), 

h4± = -2ad{W±-A). 



{D7a) 

(D76) 
(£>7c) 
{D7d) 
{Die) 



The most important factor in the above expressions is 
the first factor in the denominator of Eq. D7a, which 
can be rewritten 

W+W- - AAB = n^ujlin^ujl + 2{A + B)- A^] + {A- Bf 



Here, the first term in brackets contains the essential uj- 
dependence of the problem, and can essentially be set 
to zero in all remaining terms. 

The residual interaction terms are given by Eq. 17, 
with 



and 



k^k' ,q 



„(2) _ QqQq-gqjk) 



{D8) 



n 



-q'-q 



^ {Ak,q' P6-,k-q,q' -q ~ /2fe+g',g'^6-,ft-g,- 
'^f3k,q'Ps-,k-q,q'-q~ f 4k+q' ,q' P8-,k-q,-q' -q 
'^hk,q'Pl-,k-q,q'-q'^ hk+q',q'Pl-,k-q-q'-q 
~^'^2k—q,q'P5 — ,k^q'+q ^ Ik+q' —q^q" Pb — ^k^q—q* 
'^f 6k-q,q' Pl-,k,q'+q f 6k+q' -q,q' Pi- ,k,q-q' 
'^hk-q,q'Ps-,k,q'+q'^ hk+q'-q,q'Ps-,k,q-q')^ ^^^) 

where it is assumed that the terms linear in QqPq will 
average to zero. 



2. Canonical Transformation: Quasi-ID 
Approximation 

Due to the complicated nature of the solution, Eqs. 
Dl,7, this subsection introduces a simpler, quasi-lD ap- 
proximation: the commutator of S with the hopping 
Hamiltonian along the y-direction will be neglected. The 
next subsection gives a different, effective mass type ap- 
proximation. 

In the quasi-lD case, only the /j's and /ij's with i = 
1, 2, 5; 6 are non- vanishing. When c = d = 0, Equations 
D7 reduce to h3± = hi± = 0, /2 = 0, 



/i 



/12 = 



—ehi 
-2abhi 



and 



W+W--AAB' 



{DlOa) 
{DWb) 

(DlOc) 
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Equation DlOa can be rewritten in the suggestive form 



flk.n 



'''' ^'^}-iEn-E^,.r 



with 



and 



A /A2 



2 V 4 ^'^ 



{Dlla) 



{Dllb) 



(£>llc) 



The effective electron-electron interaction /fee can still 
be written in the form of Eq. DS, while the residual 
electron-phonon interaction becomes 



„(2) _ QqQ£gq{k)^ 



•^Ik+q^ — q,q' P5 — ,k,q— 



'^f5k,q'P2-,k-q,q'-q + f 5k+q' ,q' P2- ,k-q,-q' -q 
'^hk-q,q'Pl-,k,q'+q'^ f6k+q'-q,q'Pl-,k,q-q')- i^^'^) 

Equation D12 can be further simplified, by noting that in 

the most important case, (f = —q and 2q \s a reciprocal 
lattice vector, and hence equivalent to zero. In this case, 

rr(2) _ ^ QqQ-q9q{k) ^, 



k,q 



n 



■ {fik,-M 



qf^6-,k-q,0 ^5-,fc,0J 



'[•^6fe-g,-g~'~ J^6k,-q\Pl-,k,0 



a phase transition. A discontinuous change in pe ~ P5 
would signal an ionic-to-covalent transition. 

This canonically transformed Hamiltonian contains in- 
formation on many physical processes, not all of them rel- 
evant to the present discussion. Here, I will briefly survey 
some of the terms which arise at mean field level. Prom 

(2) 

the term in < Q^_^ > fi leads to a polaronic 
renormalization of the Cu-Ox hopping: 

(neglecting terms mk — q). The average over is given 

by 



w| < QqQ-q huJ^{N^- 



(£»156) 



where is given by Eq. 25b. There is also a renormal- 
ization of the phonon frequency, 



12 2 \ ^ 

Q Q 



V^K-M^) < P6-,k-qfi - P5-,k,0 > 

n 



(DIG) 

The presence of pe — indicates that this softening is 
associated with charge transfer effects, and will be modi- 
fied by correlation effects: in the insulating phase at half 
filling, p5 is constrained to be unity. 

In He-e, the mixed terms will lead to additional in- 
terference effects. However, the most important term is 
given by Eq. 18 in the main text. 



3. Canonical Transformation: One Band 
Approximation 

An alternative approximate solution to Eq. 11 can 
be found by assuming that only the antibonding band is 

popTilatcd by holes, and projecting the electron-phonon 
Hamiltonian onto this band. The electronic Hamiltonian 
is 



(D17) 



-[f5k,-q+ f5k-q,-q]Pl-,k-ifi)- C^^^) 

Note that the term in is of the form of Eq. 27, but 
with a quadratic Q dependence. 

In mean field theory, the four-operator terms are split 
up according to 

AB A< B > + < A> B- < Ax B >, (D14) 

where A, B are products of two operators, and < • • • > 
denotes a thermal equilibrium average. At sufficiently 
high temperatures, only QqQ-q and PjjQ, i = 5,7 have 
non-zero averages. As the temperature is lowered, other 
quantities may develop nonvanishing averages, signaling 



with Er given by Eq. 6, 



and, at the Fermi surface. 

ig{k)t{sx - Sy) 



(A)2 +4^2(^2 +,2) 



(£»186) 



In this case, by the same calculations as for the quasi- 

ID case, we find 
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^2k 



(£)20a) 



{D2Qb) 



where the primed sum means that only states which sat- 
isfy the conditions — Ep\ < Tiloq and \E^_^^ — Ep\ < 
hujo are included. In the present problem, Ep = 0. 

The resulting Hamiltonian, He + Hcdw can be diago- 
nalized by the Bogoliubov-Valatin transformation, yield- 
ing quasiparticle energies 



He-e = - E "e/2£'^K/£''^J+ (^21) 



k,k' 



(2) 

and H^_pf^ is given by Eq. 19. The phonon renormaliza- 
tion is 



^0 = ^0 - 



{D22) 



Er -\- Er - / Er — Er - 
k 2 V 2 k k 



with 



I Go otherwise. 



When E^ — E^_-., this reduces to 



4. Effective Electron-electron Interaction 



which should be compared to Eq. C8. 



(L>27c) 



The electron-electron interaction can be rewritten as 

(£)23a) 



k,k' 



with 



Pk = ^k^l^u/f^k' 



and 



Vk 



iSx {^Sx ~\~ Cx^t 



(Z3236) 



(£i23c) 



-Fo it\Ej:-Ep\<nuJo, 

and l^^fc^^j, - £^_f| < ?iwo, 
otherwise. 



{D23d) 



with Vb = AgQ/ujQ. Following Balseiro and Falicov 
this can be decoupled as 



CDW 



-GoY.Pk-G.Y.Pk + ^^ (^24) 



Ao ' 



with 



Go = Vb^ <Pfc>, 



Gi=VoY.<Pk>^ 



{D25a) 



(£1255) 



APPENDIX E: SPIN-ORBIT 
ELECTRON-PHONON COUPLING 

Earlier suggestions involving charge density waves 
(CDWs) in the cuprates were largely abandoned when 
it was pointed out that the orthorhombic splitting in 
the low-temperature orthorhombic (LTO) phase of LSCO 
leaves the two in-plane oxygens in symmetrically equiv- 
alent positions, and hence should not split the vHs de- 
generacy |7l|. In fact, there are actually two indepen- 
dent mechanisms which can lead to a splitting of the vHs 
degeneracy in the LTO phase: a dynamic Jahn- Teller 
effect, where the local symmetry is not orthorhombic, 
or spin-orbit coupling, even in a uniform orthorhombic 
phase It would be useful to be able to simplify the 
analysis by definitively ruling out one or the other effect. 
This has proven to be a non-trivial task, as the present 
estimate makes clear. 

Spin-orbit coupling can profoundly change the shape 
of the Fermi surfaces, by opening gaps where the bands 
cross, and thereby splitting the vHs degeneracy in the 
dos. This spin-orbit co upling also leads to a linear cou- 
pling to the tilt modes |72|j73[| . However, the magnitude 
of these effects is expected to be small, and hence the 
effect would be unimportant for driving either the struc- 
tural or magnetic transitions. 

Here, a simple estimate of X^p = NVq is presented, 
which reveals an underlying subtlety. The dos can be 
estimated as 



N 



16 B 

ln{ — ), 



TT^B 



(El) 



where B is the electronic bandwidth and e is an energy 
scale, ~ T. For the three-band model, B ~ 8t^/A. As- 
suming the renormalized parameters are t ~ 0.25eV, A ~ 
ley, then B ~ 0.5eV, e ~ lOmeV, and N ~ ISe^-^ 
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The electron-phonon interaction Vq can be estimated 
as follows. The interaction Hamiltonian can be written 



where the spin-orbit coupling term is \ — "f9, with 6 the 
octahedral tilt angle, 



2z 



= ri{a + a'^). 



Here, ag is the in-plane lattice constant, z the displace- 
ment of the O out of the Cu02 plane, the phonon 
creation operator, and 

mujaQ 

with m the mass of an O ion and the bare tilt-mode 
phonon frequency huj ^ bmeV. The parameter 7 is, 
approximately [7^ ] 

7 ~ — -t ~ 25meV 
9 

with Ag/g ^ 0.1 the relative g-shift. Finally, 



so 



Vn = 



A. 



AmeV. 



VqN ~ .05. 



{E2) 



This is too small to drive either the structural or su- 
perconducting transition. However, 7 ~ i, and, up to 
logarithmic corrections, TV ~ t^, so Xep ~ f^- Grilli and 
Castellani have suggested that, in the present cir- 
cumstances, Xep should not be renormalized by correlation 
effects. This would mean that the bare value t = l.SeV 
should be used in Eq. E2, leading to Xep ~ 36, a surpris- 
ingly large value. 



Table I: Parameters for Fig. 2 



Data Set io^o loq 



(1) 


11.3 


8.9 


4060 


(5) 


11.0 


8.8 


3880 


(10) 


8.8 


7.3 


2550 


(11) 


6.5 


5.9 


1500 



Fig. 2 Doping dependence of LTO transition in 
La2_2;_j^Sr2.Ndj,Cu04. Open symbols from Fig. 1 of [p7|; 
filled circles from |^^. Lines = theory: solid lines corre- 
spond to y = 0; dashed to a; = 0.15. Frames a, b, and 
c correspond to fit data sets 1, 5, and 10 respectively in 
Table I (set 11 is indistinguishible from set 10). 

Fig. 3 Phonon modes in LSCO: (a) q — qi = (7r/a,0) 
Cu-Cu bond stretching mode; (b) q = qi 0-0 bond 
stretching mode; (c) q ~ qo = (tt/a, 7r/a) Cu-Cu bond 
stretching mode; (d) q = qo 0-0 bond stretching mode; 
(e) FE mode; (f) shear strain; (g) possible mixed 0-0 
bond stretching, octahedral tilt mode. 

Fig. 4 In-plane Cu-0 bond length vs doping, after 
Radaelli, et al. |33j. Open circles = lOK data; filled cir- 
cles = 295K data; +'s = locations of the HTT-LTO tran- 
sition. Dot-dashed line = Eq. A12; other lines = guides 
to the eye. 

Fig. 5 Electronic dispersion associated with static dis- 
tortions of the symmetry of the phonon modes in Fig. 3. 
The dashed lines show the dispersion along T — Y ~ AI, 
whenever this differs from that along T — X — M, while 
the dot- (dot-) dashed lines are the dispersion along the 
line X{Y) — M = (7r/2a, 7r/2a). For ease in compari- 
son, frames e and f are shown zone-folded into the same 
V2 X supercell as in frames c,d (beaded lines show 
original bands). 

Fig. 6 Effect of correlations and polaronic effects on 
the Cu02 band dispersion in Bi-2212. Solid circles — 
data of Dessau, et al. [Q. Fig. 6a: Dotted line = bare pa- 
rameters (data for LSCO); solid line = parameters renor- 

dashed lines 



ghost 



malized by correlation effects |34 
Fermi surfaces, due to orthorhombic superlattice. Fig. 
6b: solid and dashed lines = correlated bands, with ad- 
ditional band narrowing due to polaronic effects. 

Fig. 7 Normalized Fermi velocity, v* = hirvpi/a at a 
series of Fermi levels. (Renormalized band parameters: 
t = Q.25eV, A = leV.) Fig. 7a: energies below the vHs, 
EyHs — l-207eV. Dashed lines, with increasing w*: E = 
1.01, 1.02, 1.04, 1.06, 1.08eV; solid lines, with decreasing 
V*: E = 1.10, 1.12, 1.14, 1.16, 1.18, 1.20eV, and E = 
FiyHs- Fig- 7b: energies above the vHs. Solid lines, with 
increasing w*: E = E^hs, 1-21, 1.22, 1.24, 1.26, 1.28eV; 
dashed lines, with decreasing v*: E — 1.30, 1.32, 1.33eV. 
Fig. 7c: V* (solid line) and Ej^ 
plotted vs. Ej:. 



^k-qo (dashed line) 



Figure Captions 

Fig. 1 Phase diagram of LSCO, including the pseudo- 

f transition at T = T*. Symbols from Hwang, et al. 
, short-dashed line = twice the LTO transition tem- 
perature. 



[1] D.S. Dessau, Z.-X. Shen, D.M. King, D.S. Marshall, L.W. 

Lombardo, P.H. Dickinson, A.G. Loeser, J. DiCarlo, C- 

H. Park, A. Kapitulnik, and W.E. Spicer, Phys. Rev. 

Lett. 71, 2781 (1993). 
[2] A. A. Abrikosov, J.C. Campuzano, K. Gofron, Physica 

C214, 73 (1993). 



20 



[3] J.E. Hirsch and D.J. Scalapino, Phys. Rev. Lett. 56, 2732 

(1986) ; J. Labbe and J. Bok, Europhys. Lett. 3, 1225 

(1987) ; I. E. Dzyaloshinskii, Zh. Eksp. Teor. Fiz. 93, 1487 
(1987) [Sov. Phys.-JETP 66, 848 (1987)]. For reviews, 
see J. Friedel, J. Phys. Cond. Matt. 1, 7757 (1989); R.S. 
Markicwicz, Int. J. Mod. Phys. B5, 2037 (1991); and 
D.M. Ncwns, C.C. Tsuei, P.C. Pattnaik, and C.L. Kane, 
Comments in Cond. Matt. Phys. 15. 273 (1992). 

[4] R.S. Markiewicz, Physica C193, 323 (1992). 
[5] R.S. Markiewicz, Physica C200, 65 (1992). 
[6] R.S. Markiewicz, Physica C210, 235 (1993). 
[7] R.S. Markiewicz, Physica C207, 281 (1993), and 210, 
264 (1993). 

[8] R.E. Cohen, Nature (London) 358, 136 (1992). 
[9] M. Posternak, R.Resta, and A. Baldereschi, Phys. Rev. 
B50, 8911 (1994). 
[10] T. Egami, S. Ishihara, and M. Tachiki, Science 261, 1307 

(1993) ; S. Ishihara, T. Egami, and M. Tachiki, Phys. Rev. 
B49, 8944 (1994); S. Ishihara, M. Tachiki, and T. Egami, 
Phys. Rev. B49, 16123 (1994). 

[11] W. Zhong and D. Vanderbilt, unpublished. 

[12] A.C. Hewson, "The Kondo Problem to Heavy Fermions" 
(University Press, Cambridge, 1993), p. 319; H.A. Mook, 
R.M. Nicklow, T. Ponnoy, F. Holtzbcrg, and M.W. 
Schafcr, Phys. Rev. B18, 2925 (1978); and H. Boppart, 
A. Trcindl, P. Wachter, and S. Roth, Sol. St. Commun. 
35, 483 (1980). 

[13] M. Braden, W. Reichardt, W. Schmidbauer, A.S. Ivanov, 
and A. Yu. Rumiantscv, J. Suporcond., to bo published. 

[14] F.S. Ham, in "Electron Paramagnetic Resonance", ed. 
by S. Geschwind, (N.Y., Plenum, 1972), p. 1 (esp. pp. 
52-59). 

[15] H.Y. Hwang, B. Battlogg, H. Takagi, H.L. Kao, J. Kwo, 

R.J. Cava, J.J. Krajowski, and W.F. Peck, Jr., Phys. 

Rev. Lett. 72, 2636 (1994). 
[16] R.E. Cohen, W.E. Pickett, and H. Krakauer, Phys. Rev. 

Lett. 62, 831 (1989); W.E. Pickett, R.E. Cohen, and H. 

Krakauer, Phys. Rev. Lett. 67, 228 (1991). 
[17] J. Rossat-Mignod, L.P. Rognault, P. Bourgcs, P. Burlct, 

C. Vettier, and J.Y. Henry, Physica B192, 109 (1993). 
[18] J.W. Loram, K.A. Mirza, J.R. Cooper, and W.Y. Liang, 

Phys. Rev. Lett. 71, 1740 (1993). 
[19] P. Schwciss, W. Reichardt, M. Braden, G. Collin, G. 

Hogcr, H. Glaus, and A. Erb, Phys. Rev. B49, 1387 

(1994) . 

[20] A.P. Litvinchuk, G. Thomsen, and M. Cardona, Sol. St. 

Commun. 83, 343 (1992). 
[21] X.K. Chen, E. Ahendorf, J.C. Irwin, R. Liang, and W.N. 

Hardy, Phys. Rev. B48, 10530 (1993). 
[22] H. Schwer, J. Kaxpinski, and E. Kaldis, to be published, 

Physica G (Proc. M^S-HTSC-IV). 
[23] H. Monien, D. Pines, and M. Takigawa, Phys. Rev. B43, 

258 (1991); H. Monien, P. Monthoux, and D. Pines, Phys. 

Rev. B43, 275 (1991). 
[24] B. Normand, H. Kohno, and H. Fukuyama, unpublished. 
[25] L. Tisza, personal communication. 

[26] J.B. Goodcnough and J. Zhou, Phys. Rev. B42, 4276 
(1990); J.B. Goodenough, Supercond. Sci. Technol. 3, 26 
(1990). 

[27] H. Takahashi, H. Shaked, B.A. Hunter, P.G. Radaelh, 



R.L. Hitterman, D.G. Hinks, and J.D. Jorgensen, Phys. 
Rev. B50, 3221 (1994). 

[28] B. Biichner, M. Cramm, M. Braden, W. Braunisch, O. 
Hoffels, W. Schnelle, J. Haxnischmacher, R. Borowski, 
A. Gruetz, B. Heymer, C. Hohn, R. Miiller, O. Maldon- 
ado, A. Freimuth, W. Schlabitz, G. Hogcr, D.I. Khomskii, 
and D. Wohlleben, in "Phys. and Mat. Sci. of High Tem- 
perature Superconductors, 11", cd. by R. Kossowsky, B. 
Raveau, D. Wohlleben, and S.K. Patapis, (Kluwer, Dor- 
drecht, 1992), p. 349. 

[29] J.S. Schilling and S. Klotz, in "Physical Properties of 
High Temperature Superconductors HI", ed. by D.M. 
Ginsberg (World Scientific, Singapore, 1992), p. 59, esp. 
p. 80. 

[30] J.-S. Zhou, H. Chen, and J.B. Goodenough, Phys. Rev. 

B49, 9084 (1994). 
[31] J.-S. Zhou, G.I. Bersuker, and J.B. Goodenough, to be 

published, J. Supercond. 
[32] S.J.L. Billenge and T. Egami, Phys. Rev. B47, 14386 

(1993) . 

[33] P.G. Radaelh, D.G. Hinks, A.W. Mitchell, B.L. Hunter, 
J.L. Wagner, B. Dabrowski, K.G. Vandervoort, H.K. 
Viswanathan, and J.D. Jorgensen, Phys. Rev. B49, 4163 

(1994) . 

[34] R.S. Markiewicz, Physica C228, 227 (1994). 

[35] S. Tang and J.E. Hirsch, Phys. Rev. B37, 9546 (1988). 

[36] W. Reichardt, L. Pintschovius, N. Pyka, P. Schweiss, A. 

Erb, P. Bourges, G. Collin, J. Rossat-Mignod, I.Y. Henry, 

A.S. Ivanov, N.L. Mitrofanov, and A.Yu. Rumiantsev, J. 

Supercond. 7, 399 (1994). 
[37] G. Bilbro and W.L. McMillan, Phys. Rev. B14, 1887 

(1976). 

[38] A.M. Gabovich, D.P. Moiseev, and A.S. Shpigel, J. Phys. 

C15, L569 (1982). 
[39] J. Song and J.F. Annett, unpublished. 
[40] P.C. Hammel, A.P. Reyes, S-W. Cheong, Z. Fisk, and 

J.E. Schirber, Phys. Rev. Lett. 71, 440 (1993). 
[41] S.J.L. Billinge, G.H. Kwei, and H. Takagi, Phys. Rev. 

Lett. 72, 2282 (1994). 
[42] D.J. Scalapino, R.T. Scalettar, and N.E. Bickers, in 

"Novel Superconductivity", ed. by S.A. Wolf and V.Z. 

Krcsin, (Plenum, NY, 1987), p. 475. 
[43] A. A. Aligia, M. Kulic, V. Zlatic, and K.H. Bcrmcmann, 

Sol. St. Commun. 65, 501 (1988). 
[44] J.H. Kim and Z. Tesanovic, Phys. Rev. Lett. 71, 4218 

(1994). 

[45] R.S. Markicwicz, unpublished. 
[46] E. Pytte, Phys. Rev. BIO, 4637 (1974). 
[47] C.A. Balseiro and L.M. FaUcov, Phys. Rev. B20, 4457 
(1979). 

[48] J. Bardeen and D. Pines, Phys. Rev. 99, 1140 (1955). 

[49] T. Holstcin, Ann. Phys. 8, 325, 343 (1959). 

[50] G.D. Mahan, "Many-Particle Physics" (N.Y., Plermm, 

1981), Sect. 6-2. 
[51] B.I. Halperin and T.M. Rice, in "Solid State Physics", 

Vol. 21, Ed. by F. Seitz, D. TurnbuU, and H. Ehrenreich 

(NY, Academic, 1968), p. 115. 
[52] H.J. Schulz, Phys. Rev. B39, 2940 (1989), and A.A. Ners- 

esyan, G.I. Japaridze, and I.G. Kimeridze, J. Phys. Cond. 

Matt. 3, 3353 (1991). 



21 



[53] I. Affleck and J.B. Marston, Phys. Rev. B37, 3774 (1988), 

G. Kotliar, Phys. Rev. B37, 3664 (1988). 
[54] M. Grilli and C. Castellani, Phys. Rev. B50, 16880 

(1994). 

[55] A.C. Hewson and D.M. Newns, J. Phys. C 12, 1665 

(1979). 

[56] D. Sherrington and P. Riseborough, J. de Phys. C-4, 255 
(1976). 

[57] B. Alascio and A. Lopez, Sol. St. Commun. 14, 321 
(1974). 

[58] P.W. Anderson and S.T. Chui, Phys. Rev. B9, 3229 
(1974). 

[59] A.S. Alexandrov, J. Ranninger, and S. Robaszkiewicz, 

Phys. Rev. B 33, 4526 (1986); A.S. Alexandrov and N.F. 

Mott, J. Supercond. 7, 599 (1994); D. Emin and M.S. 

Hillery, Phys. Rev. B39, 6575 (1989). 
[60] R.D. Shannon, Acta Cryst. A32, 751 (1976). 
[61] I.D. Brown and D. Altermatt, Acta Cryst. B41, 244 

(1985). 

[62] R.M. Hazeri, in "Physical Properties of High- Tempera- 
ture Superconductors 11", ed. by D.M. Ginsberg (Singa- 
pore, World, 1990), p. 121. 

[63] L. Pauling, J. Am. Chem. Soc. 49, 765 (1927). 

[64] M. Brcuor, B. Biichncr, R. Miillor, M. Crarnrn, O. 
Maldonado, A. Freimuth, B. Roden, R. Borowski, B. 
Heymer, and D. Wohlleben, Physica C208, 217 (1993). 

[65] W. Ting, K. Fossheim, and T. Laegreid, Sol. St. Com- 
mun. 80, 47 (1991). 

[66] D.M. dcLccuw, W.A. Grocn, L.F. Feiner, and E.E. 
Haviuga, Physica C166, 133 (1990). 

[67] ,J.L. Tallon and G.V.M. WilUams, J. Less-Comm. Mets. 
164-5, 70 (1990). 

[68] S. Tanaka, N. Fukushima, H. Niu, and K. Ando, Jpn. J. 
Appl. Phys. 29, L1987 (1990). 

[69] J.C. Grenier, presented at M^S-HTSC-IV, Grenoble, 
France, July 5-9, 1994, and personal communication. 

[70] H. Bilz, G. Benedek, and A. Bussmann-Holder, Phys. 
Rev. B35, 4840 (1987); A. Bussmann-Holder, A. Simon, 
and H. Biittncr, Phys. Rev. B39, 207 (1989). 

[71] J.P. Pouget, C. Noguera, and R. Moret, J. Phys. (France) 
49, 375, 1988. 

[72] D. Coffey, T.M. Rice, and E.G. Zhang, Phys. Rev. B44, 

10112 (1991). 

[73] N.E. Bonesteel, T.M. Rice, and E.G. Zhang, Phys. Rev. 
Lett. 68, 2684 (1992). 



22 



